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A  physical  and  numerical  model  of  two-stage  thermoelectric  energy  harvesting  system  driven  by  blast 
furnace  slag  water  waste  heat  is  established.  The  performance  of  the  system  with  counter-flow  type  heat 
exchangers  is  investigated  by  numerical  simulation.  In  the  case  of  the  temperatures  of  heat  reservoirs 
change  over  flow  passage,  the  effects  of  inlet  temperature  of  flushing  slag  water,  convective  heat  transfer 
coefficient  and  flow  passage  length  on  the  power  output,  efficiency,  maximum  power  output  and 
maximum  efficiency  as  well  as  optimal  resistance  ratio  of  the  system  are  analyzed.  Moreover,  the 
electrical  current  range  corresponding  to  the  maximum  power  output  and  maximum  efficiency  is  ob¬ 
tained.  Simulation  results  show  that  the  maximum  power  output  of  0.44  kW  and  maximum  efficiency  of 
2.66%  are  available  with  inlet  temperature  of  blast  furnace  slag  water  at  100  °C  if  load  resistance  is 
matched.  The  optimal  resistance  ratio  corresponding  to  the  maximum  power  output  is  about  1.13. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Thermoelectric  generation  is  classified  as  direct  power  conver¬ 
sion  [1—5],  It  has  lots  of  features  and  advantages  compared  with 
traditional  generators,  for  example,  the  absence  of  moving  com¬ 
ponents  results  in  an  increase  of  reliability,  and  reduction  of  extra 
maintenance;  the  modularity  can  provide  a  wide-scale  application 
without  significant  losses  in  performance.  Besides,  these  devices 
produce  no  noise  and  waste  in  the  conversion  process.  Therefore, 
thermoelectric  generator  has  been  regarded  as  a  useful  and 
attractive  device  for  direct  energy  conversion  [6-10]  and  potential 
application  in  waste  heat  recovery  [11-14], 

As  energy  shortage  and  environmental  deterioration  are 
growing,  energy  saving  and  emission  reduction  have  become  global 
obligation  and  responsibility  [15-19],  Yu  and  Zhao  [15]  employed 
hot  and  cold  water  as  heat  reservoirs  and  analyzed  the  performance 
of  thermoelectric  generator  with  the  parallel-plate  heat  exchanger. 
Meng  et  al.  [16]  established  a  numerical  model  of  commercial 
thermoelectric  generator  with  finned  heat  exchangers  taking  into 
account  inner  and  external  multi-irreversibilities.  Gou  et  al.  [17] 
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established  a  dynamic  model  for  waste  heat  recovery  in  thermo¬ 
electric  generator  to  assess  the  effects  of  heat  reservoirs  on  the 
dynamic  characteristics.  Kazuaki  et  al.  [18]  analyzed  energy  econ¬ 
omy  for  a  combined  thermoelectric  generator  on  top  of  a  steam 
turbine  cycle,  and  demonstrated  the  advantage  of  adding  a  ther¬ 
moelectric  on  top  of  a  steam  turbine  cycle.  Stevens  [  19]  established 
an  energy  harvesting  device  that  produces  milliwatt-scale  power 
uses  a  thermoelectric  generator  operating  between  the  air  and 
ground  temperatures. 

Since  the  single  stage  thermoelectric  generator  cannot  operate 
in  the  case  of  large  temperature  range  [20-25],  and  due  to  the 
performance  limits  of  thermoelectric  materials  and  requirements 
of  adapting  different  heat  reservoir  conditions,  some  authors  have 
investigated  the  performance  of  thermoelectric  generator 
composed  of  two  stages  and  more  [26-32], 

In  the  process  of  promoting  the  transformation  of  materials  in 
the  iron  and  steel  industry,  rich  and  variety  of  waste  heat  is 
generated.  Taking  advantage  of  these  waste  heats  can  effectively 
reduce  energy  consumption  [33—35],  Chen  et  al.  [36]  provided  an 
overview  of  residual  heat  recovery  in  iron  and  steel  enterprises  in 
China  based  on  single  stage  and  multi-element  thermoelectric 
generation,  and  the  calculations  showed  that  electricity  power  of 
21  kWh,  43  kWh  and  60  kWh  can  be  recovered  from  1  GJ  waste 
heat  with  temperature  difference  of  100  °C,  200  °C  and  300  °C, 
respectively.  Meng  et  al.  [37]  established  a  model  of  single  stage 
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Nomenclature 

Greek  symbol 

a 

Seebeck  coefficient  (V  K_1) 

A 

area  (m2) 

a 

electrical  conductivity  (Q_1  m-1) 

B 

flow  channel  width  (m) 

A 

thermal  conductivity  (W  m-1  K-1) 

c 

specific  heat  capacity  at  constant  pressure  (J  kg-1 1<  1 ) 

e 

fill  factor 

c 

mass  flow  rate  (kg  s-1) 

V 

conversion  efficiency 

H 

flow  channel  height  (m) 

h 

convective  heat  transfer  coefficient  (W  m-2  K-1) 

Subscripts 

I 

electrical  current  (A) 

c 

cold  junction 

K 

thermal  conductance  (W  K-1) 

f 

flow  passage 

L 

length  (m) 

g 

air  gap 

P 

power  output  (W) 

h 

hot  junction 

a 

heat  flow  rate  (W) 

L 

load 

R 

electrical  resistance  (Q) 

m 

intermediate  junction 

r 

load  resistance  ratio 

1 

heat  source 

T 

thermodynamic  temperature  (K) 

2 

heat  sink 

1 

Celsius  temperature  (°C) 

and  multi-element  thermoelectric  generator  driven  by  blast 
furnace  slag  flushing  water  sensible  heat  with  parallel-flow  type 
heat  exchangers,  and  analyzed  the  power  and  efficiency  perfor¬ 
mance,  as  well  as  the  recovery  period  of  the  equipment  cost. 

On  the  basis  of  research  achievements  mentioned  above,  this 
paper  will  establish  a  physical  and  numerical  model  of  two  stage 
thermoelectric  energy  harvesting  system  driven  by  blast  furnace 
slag  water  waste  heat,  and  analyze  the  effects  of  key  parameters  on 
the  power  output,  efficiency,  optimal  electrical  current  as  well  as 
optimal  resistance  ratio  of  the  system. 


2.  Harvesting  system  setup  and  model 

A  schematic  diagram  of  the  harvesting  system  driven  by  blast 
furnace  slag  water  waste  heat  is  shown  in  Fig.  1.  The  system  is 
consisted  of  two  parts.  The  first  part  contains  two  heat  exchangers 
between  heat  reservoir  and  two-stage  thermoelectric  generator 
module.  As  is  well  known,  the  performance  of  counter-flow  heat 
exchanger  is  better  than  that  of  parallel-flow  one,  this  paper  aims  at 
discussing  the  system’s  performance  with  counter-flow  type  heat 
exchangers  between  blast  furnace  slag  water  and  cooling  water. 
The  second  part  is  the  two-stage  thermoelectric  generator.  The 
waste  slag  washing  water  flows  through  the  hot-side  heat- 


exchanger  channels  and  heats  the  upper  surface  of  thermoelectric 
module.  The  cooling  water  driven  by  a  pump  flows  through  the 
cold-side  heat-exchanger  channels  and  absorbs  heat  from  the 
lower  surface  of  thermoelectric  modules. 

The  generator  consisted  of  a  top  stage  and  a  bottom  stage  with 
the  same  pairs  of  thermoelectric  elements.  Each  element  is 
composed  of  p-type  and  n-type  semiconductor  legs.  The  thermo¬ 
electric  power-generation  element  is  assumed  to  be  insulted,  both 
electrically  and  thermally,  from  its  surroundings,  except  at  the 
junction-reservoir  contacts  and  the  junction  between  the  two 
stages.  The  structure  parameters  of  two  heat-exchangers  and  the 
geometric  dimensions  of  each  element  of  the  two-stage  modules 
are  assumed  to  be  the  same. 

The  thermodynamic  model  of  two-stage  thermoelectric  gener¬ 
ator  system  driven  by  blast  furnace  slag  waste  heat  is  shown  in 
Fig.  2.  The  temperatures  of  slag  washing  water  and  cooling  water 
are  Ti(x)  and  T2(x),  respectively.  The  temperatures  of  the  hot- 
junction  and  cold-junction  are  Th(x)  and  Tc(x),  respectively.  The 
top  stage  and  bottom  stage  of  the  thermoelectric  generator  are 
assumed  to  be  well  contacted,  so  the  thermal  contact  resistance 
between  the  two  modules  is  neglected.  The  temperature  of  the 
cold-side  junction  of  top  stage  and  the  temperature  of  the  hot-side 
junction  of  bottom  stage  are  assumed  to  be  the  same,  i.e.  Tm(x). 

Since  the  thermoelectric  thermocouples  are  not  closely  packed 
for  insulation,  there  exists  an  air  gap  between  the  thermocouples. 
Thus  heat  flow  from  the  heat  reservoir  directly  flows  through  the 
air  gap  rather  than  through  the  thermocouples,  and  the  heat  flow  is 


T/x) 


waste  heat. 
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expressed  by  Qg.  Taking  into  account  the  heat  flow  through  the  air 
gap  between  the  thermocouples,  the  fill  factor  6  (0  <  d  <  1)  [22]  is 
introduced  to  express  the  proportion  of  the  thermoelectric  ele¬ 
ments  area  on  the  modules: 


8  =  2AN/Acp  (1) 

where  A,  N  and  Acp  are,  respectively,  the  cross-sectional  area  of 
thermo-elements,  the  number  of  thermocouples  in  a  module,  and 
the  cross-sectional  area  of  a  module.  The  heat  flow  rate  of  the  air 
gap  leakage  of  top  stage  generator  and  bottom  stage  generator  is 
Qgi  and  Qg2,  respectively. 

For  the  top  stage  generator,  the  heat  flow  rate  absorbed  from  the 
heat  source  is  Qi.  For  the  bottom  stage  generator,  the  heat  flow  rate 
dissipated  to  the  heat  sink  is  Q2.  Moreover,  the  rates  of  heat  ab¬ 
sorption  and  dissipation  at  hot  junction  of  top  stage  generator  are 
Qhi  and  Qci,  respectively.  The  rates  of  heat  absorption  and  dissi¬ 
pation  at  cold  junction  of  bottom  stage  generator  are  Qh2  and  Qc2. 
respectively.  The  two  stage  thermoelectric  generator  is  connected 
electrically  in  series  along  the  directions  of  fluid  path,  but  in  parallel 
for  the  directions  perpendicular  to  fluid  path.  The  power  output  of 
the  system  is  P.  The  electrical  resistance  of  the  external  load  and 
single  element  are  Ri  and  R,  respectively.  The  electrical  current 
generating  through  the  system  at  steady  condition  is  I.  The  heat 
conductances  of  hot  and  cold  side  heat-exchanger  are  K\  and  K2, 
and  both  of  them  are  calculated  from  four  thermal  resistances:  heat 
exchanger  separator  resistances  oex//ex,  ceramic  substrate  re¬ 
sistances  (5Cp/ACp  and  convective  heat  transfer  resistances  1  //icv  as 
well  as  heat  conduction  grease  contact  resistances  5C/A c-  That  is: 


Kl  ~  WAex  +  acp/^  +  l/ftcv  +  ac/V  Kl  “ /Cl  (2) 

A  control  volume  shown  in  Fig.  2  is  formulated  for  a  set  of  ele¬ 
ments,  which  includes  flow  passages  on  both  sides  of  the  ther¬ 
moelectric  generator.  Thus,  the  governing  equations  of  the  control 
volume  for  the  top  stage  generator,  respectively,  are 


dQi  =  -GiCpi  d7j  (x) 


dQi  =K,[Mx)-Th(x)]B  dx 


dQhi  =  [«rh(x)/  +  K(Tb(x)  -  Tm(x))  - 


dQci  =  [aTm(x)f +  K(rh(x)  -  Tm(x))  + 


(3) 

(4) 

I2R1 

6B  dx 

(5) 

2  J 

2A 

I2R] 

6B  dx 

(6) 

tI 

2A 

[J2D1  0D  /jy 

aTc(x)/  +  /C(Tm(x)-Tc(x))+^-J  (10) 

where  G2,  cP2  are,  respectively,  the  rate  of  cooling  water  flow, 
specific  heat  at  constant  pressure  of  cooling  water. 

The  heat  flow  rates  of  the  air  gap  leakage  of  top  stage  generator 
and  bottom  stage  generator  are,  respectively: 

dQgi  =  Ksi  (Th(x)  -  rm(x))(l  -  6)B  dx  (11) 

d0g2  =  Kg2(rm(x)  -  TC(J0)(1  -  0)B  dx  (12) 

where  K%  1  and  K% 2  are,  respectively,  the  heat  conductance  of  the  air 
gap  leakage  of  top  stage  generator  and  bottom  stage  generator. 
Kgi  =  Kg 2  —  Aair/T,  where  Aair  is  the  thermal  conductivity  of  air,  and  L 
is  the  length  of  a  single  element. 

Based  on  the  energy  balance,  one  has: 

Qi  =  Qhi  +  Qgi  (13) 

Qm  =  Qcl  +  Qgi  =  0h2  +  0g2  (14) 


02  =  0c2  +  0g2 


Substituting  Eqs.  (3) — (12)  into  Eqs.  (13)— (15)  yields  the  differential 
equations  about  Ti(x),  Th(x),  Tm(x),  Tc(x)  and  T2(x)  as 


—GiCpi 


+  Kgi  (Th(x)  -  Tm(x))(l  -  6) 

[Kgi  (Mx)  -  rm(x))  -  %(Tm(x)  -  rc(x))]  (1  -  6) 
W  [K(2Tm(x)  -  Tc(x)  -  Th(x))  -  I2R]  A 


K2(Tc(x)  -  T2(x))  =  [«Tc(x)/  +  K(Tm(x)  -  Tc(x))  ±  ^ 


Mx)  = 


+  Kg2(Tm(x)-Tc(X))(t- 

K2(Tc(x)  -  T2(x)]B 
-C2cp2 


where  Gi,  cpi,  B,  a  and  I<  are,  respectively,  the  rate  of  blast  furnace 
slag  water  flow,  specific  heat  at  constant  pressure  of  blast  furnace 
slag  water,  width  of  fluid  passage,  Seebeck  coefficient,  thermal 
conductance. 

Making  the  same  analyses  for  bottom  stage  generator,  one  has: 


dQ2  =  -G2cp2  dT2(x) 

(7) 

dQ2  =  K2[Tc(x)  -  T2(x)]B  dx 

(8) 

The  boundary  conditions  of  temperatures  of  fluids  for  counter¬ 
flow  type  heat  exchangers  are: 

T,(0)=T10  (21) 

T2(Lf)=T20  (22) 

where  Lf  is  the  length  of  the  flow  passages  for  two  heat-exchangers. 

The  power  output  P  of  the  system  is  the  difference  between  the 
absorbed  and  dissipated  heat  flow  rate,  and  the  thermal  efficiency 
of  the  system  is  77  =  P/Qp  Then,  one  has 


(9)  P  =  G1Cp1[r1(0)-r1(Lf)] -C2Cp2[T2(0)-r2(Lf)]  (23) 


d0h2  =  [aTm(x)I  +  K(Tm(x)  -  Te(x))  ^ 
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_<*.[r,(») -*(.,)]  (24, 

Cicpl  [T,(0)-T,(Lf)\ 

The  differential  equations  (16)— (20)  belong  to  the  two-point 
boundary  value  problems,  which  cannot  be  directly  solved  with  or¬ 
dinary  methods.  In  this  paper,  trial  method  is  used  to  translate  the 
boundary  value  problem  into  initial  value  problem  that  can  be  solved 
by  numerical  method.  The  calculation  flow  chart  of  differential 
equations  is  shown  in  Fig.  3.  The  initial  T2(0)  is  input  for  calculating 
the  differential  equations,  and  the  accurate  T2(0)  can  be  determined 
through  the  judge  and  loop  statements.  Thus,  the  temperatures  that 
vary  with  the  flow  passage,  i.e.  Ti(x),  Th(x),  Tm(x),  Tc(x)  and  T2(x)  can  be 
calculated,  the  analytical  formula  describing  the  power  output  and 
efficiency  versus  the  key  parameters  are  derived. 

3.  Numerical  examples  and  analyses 

The  effects  of  some  key  parameters  on  the  system's  performance 
are  investigated  by  numerical  simulation.  Since  the  inlet  temper¬ 
ature  of  blast  furnace  slag  water  at  constant  pressure  is  less  than 
100  °C  and  the  inlet  temperature  of  cooling  water  is  set  as  20  °C,  the 
effect  of  temperature  on  the  physical  properties  of  the  thermo¬ 
electric  materials  is  ignored.  The  physical  properties  of  the  ther¬ 
moelectric  materials  at  mean  temperature  of  50  °C  are  set  as 
follows  [38]:  a  =  2.07  x  10"4  V  K"1,  A  =  1.53  W  m”1  K-1,  and 
a  =  8.70  x  104  m  Q-1.  The  length  and  cross-sectional  area  of  a 
thermoelectric  semiconductor  leg  are  L  =  2.0  x  10  3  m  and 
A=1.2x  1.2  x  10-6  m2,  respectively,  the  fill  factor  can  be  calculated 
as  6  =  0.5.  The  plate  heat  exchangers  are  made  of  aluminum.  The 
thickness  of  heat  exchangers  is  5ex  =  2.0  x  10-3  m,  the  width  of  the 
flow  channel  is  B  =  0.1  m,  and  the  height  of  the  flow  channel  is 
H=  6.0  x  10-3  m.  The  specific  heats  of  blast  furnace  slag  water  and 
cooling  water  at  constant  pressure  are  cpi  =  4.2  x  103  J  kg-1  K-1  and 
cP2  =  4.18  x  103  J  kg  1  K-1,  respectively.  The  mass  flow  rates  and 
flow  velocity  of  water  are  set  as  Gi  =  C2  =  0.6  kg  s_1  and 
V!  =  v2  =  1.0  m/s,  respectively.  Besides,  since  the  estimated 
pumping  power  (several  watts)  is  quite  small  compared  with  the 
power  output  of  the  systems,  the  effect  of  pump  power  and  its 
energy  consumption  is  not  taken  into  account. 


Fig.  3.  Calculation  flow  chart  of  differential  equations. 


Fig.  4.  Effect  of  inlet  temperature  of  flushing  slag  water  on  power  versus  electrical 
current 


3.1.  The  effects  of  inlet  temperature  of  blast  furnace  slag  water 

Figs.  4  and  5  show  the  effects  of  blast  furnace  slag  water  inlet 
temperature  on  power  output  and  efficiency  versus  working  elec¬ 
trical  current,  respectively.  Since  the  two  heat  reservoirs  are  driven 
by  water  pump  and  defined  as  forced  convective  heat  transfer,  the 
convective  heat  transfer  coefficient  of  water  is  calculated  as 
hcv  =  3000  W  m-2  K-1  based  on  the  heat  exchanger  geometry.  It  can 
be  seen  that  the  power  output  and  efficiency  versus  working  elec¬ 
trical  current  characteristic  behaves  a  parabola-like  for  fixed  inlet 
temperature  of  blast  furnace  slag  water.  There  are  two  optimal 
working  electrical  currents  Ip  and  Iv  corresponding  to  the  maximum 
power  output  and  maximum  efficiency,  respectively.  Both  the  power 
output  and  efficiency  increase  with  the  increase  of  inlet  temperature 
of  blast  furnace  slag  water.  The  maximum  power  output  of  0.44  kW 
and  maximum  efficiency  of  2.66%  are  available  for  inlet  temperature 
of  blast  furnace  slag  water  at  100  °C  if  load  resistance  is  matched. 

Fig.  6  shows  the  variation  of  maximum  power  output  and 
maximum  efficiency  with  the  inlet  temperature  of  blast  furnace 
slag  water.  It  can  be  seen  that  both  the  maximum  power  output  and 
maximum  efficiency  increase  with  the  inlet  temperature  of  blast 
furnace  slag  water.  Moreover,  the  maximum  efficiency  is  a  linear 
increasing  function  of  the  inlet  temperature  of  blast  furnace  slag 
water.  However,  the  maximum  power  output  is  a  nonlinearly 
increasing  function  of  inlet  temperature  of  blast  furnace  slag  water 
and  the  growth  rate  is  larger  and  larger.  Fig.  7  shows  the  variation  of 


Fig.  5.  Effect  of  inlet  temperature 
current. 
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Fig.  6.  The  maximum  power  and  maximum  efficiency  versus  inlet  temperature  of 
flushing  slag  water. 

the  optimal  working  electrical  currents  with  the  inlet  temperature 
of  blast  furnace  slag  water.  It  can  be  seen  that  both  /P  and  Iv  are 
linear  increasing  functions  of  the  inlet  temperature  of  blast  furnace 
slag  water,  and  the  optimal  working  electrical  current  (fp)  corre¬ 
sponding  to  the  maximum  power  output  is  larger  than  that  (/,,) 
corresponding  to  the  maximum  efficiency  for  fixed  inlet  tempera¬ 
tures.  The  intermediate  range  between  two  lines  is  the  optimal 
working  region  of  electrical  current.  Therefore,  an  appropriate 
working  electrical  current  should  be  determined  in  order  to  meet 
different  needs  in  practical  applications. 

The  inlet  temperature  of  blast  furnace  slag  water  as  industrial 
circulating  water  at  constant  pressure  is  less  than  100  °C.  However, 
it  is  suggested  that  hot  water/glycol  mixture  fluids  could  be  an 
appropriate  choice  for  the  temperature  over  100  °C  in  practical 
application  [15].  Calculation  results  show  that  the  maximum  power 
output  and  maximum  efficiency  of  the  system  can,  respectively, 
reach  to  1.14  kW  and  4.20%  for  inlet  temperature  of  blast  furnace 
slag  water  at  150  °C,  they  are  respectively  162%  and  57.8%  im¬ 
provements  than  those  of  the  system  with  inlet  temperature  of 
blast  furnace  slag  water  at  100  °C.  To  sum  up,  the  increase  of  inlet 
temperature  of  blast  furnace  slag  water  can  significantly  improve 
the  performance  of  the  system. 

3.2.  The  effects  of  convective  heat  transfer  coefficient 

Figs.  8  and  9  show  the  effects  of  convective  heat  transfer  coef¬ 
ficient  on  power  output  and  efficiency  versus  working  electrical 
current,  respectively.  Fig.  10  shows  the  variation  of  maximum  po¬ 
wer  output  and  maximum  efficiency  versus  the  convective  heat 


Working  electrical  current  (A) 

Fig.  8.  Effect  of  convective  heat  transfer  coefficient  on  power  versus  electrical  current. 

transfer  coefficient.  In  the  calculations,  Lf  =  10  m  and  t10  -  100  °C 
are  set.  It  can  be  seen  that  both  the  power  output  and  efficiency 
increase  with  the  increase  of  the  convective  heat  transfer  coeffi¬ 
cient.  However,  the  growth  rate  is  smaller  and  smaller  as  the  in¬ 
crease  of  the  convective  heat  transfer  coefficient.  The  convective 
heat  transfer  coefficient  of  water  is  about  3000-5000  W  m-2  K-1, 
the  maximum  power  output  of  0.49  kW  and  maximum  efficiency  of 
2.81%  are  available  for  inlet  temperature  of  blast  furnace  slag  water 
at  100  °C.  When  the  external  heat  resistance  losses  between  the 
thermoelectric  generator  and  the  external  reservoir  are  ignored,  i.e. 
hcv  -►  oo,  the  maximum  power  output  and  maximum  efficiency  are 
0.58  kW  and  3.06%,  respectively. 

Figs.  11  and  12  show  the  optimal  working  electrical  current  and 
optimal  resistance  ratio  versus  convective  heat  transfer  coefficient, 
respectively.  It  can  be  seen  that  the  optimal  working  electrical 
current  /opt  increases,  while  the  optimal  resistance  ratio  ropt  de¬ 
creases  with  the  increase  of  convective  heat  transfer  coefficient  at 
the  beginning.  However,  both  lopt  and  ropt  tend  to  be  stable  as  the 
convective  heat  transfer  coefficient  increases.  Moreover,  /P  >  Iv  and 
rv  >  rP  are  guaranteed  for  a  fixed  convective  heat  transfer  coeffi¬ 
cient.  The  optimal  resistance  ratio  corresponding  to  the  maximum 
power  output  is  about  1.13  and  that  corresponding  to  the  maximum 
efficiency  is  about  1.42  with  the  convective  heat  transfer  coefficient 
of  3000  W  m-2  K-1. 

It  should  be  pointed  out  that  most  of  previous  researches  about 
the  thermoelectric  generator  were  performed  by  using  non¬ 
equilibrium  thermodynamics  without  considering  the  external 


Fig.  9.  Effect  of  convective  heat  transfer  coefficient  on  efficiency  versus  electrical 
current 
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coefficient. 


Convective  heat  transfer  coefficienct  (103Wm"2K“1) 

Fig.  12.  Optimal  resistance  ratio  versus  convective  heat  transfer  coefficient. 


heat  resistance  losses  and  were  based  on  the  infinite  thermal  ca¬ 
pacity,  i.e.  constant  temperature  heat  reservoirs.  In  this  case,  the 
load  resistance  is  selected  to  equal  to  the  effective  internal  resis¬ 
tance  of  the  thermoelectric  generator  so  that  the  maximum  power 
output  of  the  system  could  be  achieved  [15,37],  i.e.  rp  =  1.  However, 
the  heat  reservoirs  in  practical  application  are  always  finite- 
capacity  and  variable-temperature.  Thus,  the  temperatures  of  the 
modules  vary  with  some  key  factors,  such  as  working  electrical 
current  and  convective  heat  transfer  coefficient.  Calculation  results 
show  that  the  optimal  resistance  ratio  corresponding  to  the 
maximum  power  output  is  rp  >  1.  Therefore,  the  external  heat 
resistance  losses  cannot  be  ignored. 

3.3.  The  effects  of  flow  passage  length 

Fig.  13  shows  the  effects  of  flow  passage  length  on  the  maximum 
power  output  and  maximum  efficiency,  respectively.  Fig.  14  is 
employed  to  show  how  the  power  per  area  changes  with  the  flow 
passage.  In  the  calculations,  tio  =  100  °C  and  ha  —  3000  W  m-2  K-1 
are  set.  It  can  be  seen  that  the  maximum  power  output  increases 
with  the  increase  of  flow  passage  length,  while  both  the  maximum 
efficiency  and  maximum  power-per-area  are  linear  decreasing 
functions  of  the  flow  passage  length.  It  results  from  the  decrease  of 
temperature  of  blast  furnace  slag  water  along  the  flow  passage  and 
then  the  temperature  difference  through  the  modules. 

Fig.  15  shows  the  effect  of  flow  passage  length  on  the  optimal 
working  electrical  current.  It  can  be  seen  that  the  two  optimal 
working  electrical  currents  corresponding  to  the  maximum  power 


output  and  the  maximum  efficiency  are  linear  decreasing  functions 
of  the  flow  passage  length.  Moreover,  lP  >  Iv  is  guaranteed  for  fixed 
flow  passage  length.  The  intermediate  range  between  two  lines  is 
the  optimal  working  region  of  electrical  current. 

4.  Conclusion 

This  paper  establishes  a  model  of  two-stage  thermoelectric 
energy  harvesting  system.  The  performance  of  the  two-stage 
thermoelectric  generator  system  driven  by  blast  furnace  slag  wa¬ 
ter  waste  heat  is  analyzed.  The  major  results  are  as  follows: 

(1)  The  inlet  temperature  of  flushing  slag  water  is  an  important 
parameter  which  affects  the  performance  of  the  two-stage 
thermoelectric  generator  system  significantly.  Both  the 
maximum  efficiency  and  the  optimal  working  electrical 
current  are  linear  increasing  functions  of  the  flushing  slag 
water  inlet  temperature,  while  the  maximum  power  output 
is  nonlinear. 

(2)  The  optimal  resistance  ratio  corresponding  to  the  maximum 
power  output  is  about  1.13  and  the  optimal  resistance  ratio 
corresponding  to  the  maximum  efficiency  is  1.42  in  the  case 
of  external  heat  resistance  losses  between  the  device  and 
heat  reservoirs  are  taking  into  account. 

(3)  Both  the  maximum  efficiency  and  the  growth  rate  of  the 
maximum  power  output  decrease  with  increase  of  flow 
passage  while  the  maximum  power  output  increases. 
Moreover,  the  two  optimal  working  electrical  currents 
decrease. 


versus  convective  heat  transfer  coefficient. 


Fig.  13.  The  maximum  power  and  maximum  efficiency  versus 


Fig.  11.  Optimal  electrical  currents ' 


flow  passage  length. 
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